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We calculate the inclusive and diffractive photoproduction of heavy quarks in proton-proton 
collisions at Tevatron and LHC energies, where the photon reaches energies larger than those ones 
accessible at DESY-HERA. The integrated cross section and the rapidity distributions for charm and 
bottom production are computed within the color dipole picture employing three phenomenological 
saturation models based on the Color Glass Condensate formalism. Our results demonstrate that 
the experimental analyzes of these reactions is feasible and that the cross sections are sensitive to 
the underlying parton dynamics. 

PACS numbers: 25.75.Dw, 13.60.Lc 

I. INTRODUCTION 

The cross sections for heavy quark production in hadron-hadron and lepton-hadron collisions at high energies are 
strongly dependent on the behavior of the gluon distribution, which is determined by the underlying QCD dynamics 
(See e.g. Refs. [H, [1]). Theoretically, at high energies the QCD evolution leads to a system with high gluon density, 
characterized by the limitation on the maximum phase-space parton density that can be reached in the hadron 
wavefunction (parton saturation). The transition is specified by a typical scale, which is energy dependent and is 
called saturation scale Q sa t (For recent reviews see Ref. Q). Signals of parton saturation have already been observed 
both in ep deep inelastic scattering at HERA and in deuteron-gold collisions at RHIC (See, e.g. Ref. j4]). In 
particular, in Ref. [|| we demonstrated that the inclusive charm total cross section exhibits the property of geometric 
scaling , which is one of the main characteristics of the high density approaches. However, the observation of this 
new regime still needs confirmation and so there is an active search for new experimental signatures. In this paper 
we study the inclusive and diffractive photoproduction of heavy quarks in proton-proton collisions considering three 
phenomenological models based on the Color Glass Condensate, which describe quite well the current experimental 
HERA data for inclusive and exclusive observables. Our goal is twofold: update our previous studies [1,0] considering 
these new parameterizations for the dipole scattering amplitude and present a comparison between the inclusive and 
diffractive production mechanisms using an identical theoretical input. 

Our main motivation comes from the fact that in coherent interactions at Tevatron and LHC the photon reaches 
energies higher than those currently accessible at DESY - HERA. In hadron-hadron colliders, the relativistic protons 
give rise to strong electromagnetic fields, which can interact with each other. Namely, quasi-real photons scatters off 
protons at very high energies in the current hadron colliders (For recent reviews on coherent interactions see, e.g., 
Ref. [|[). In particular, the heavy quark photoproduction cross section in a proton-proton collision is given by, 

a(p + p^p + QQ + Y) = 2 f ^^a lp ^ Q Q Y (w* p = 2w^S^) du , (1) 



where uj is the photon energy in the center-of-mass frame (c.m.s.), W lp is the c.m.s. photon-proton energy and \/Snn 
denotes the proton-proton c.m.s. energy. The final state Y can be a hadronic state generated by the fragmentation of 
one of the colliding protons (inclusive production) or a proton (diffractive production) . The photon spectrum is given 
by!, 

with the notation il = 1 + [ (0.71 GeV 2 )/<2 2 nin ] and Q 2 min = w 2 /[t£ (1 - 2u/VSnn)} ~ (^/7l) 2 , where j L is the 
Lorentz factor. The expression above is derived considering the Weizsacker- Williams method of virtual photons and 
using an elastic proton form factor (For more detail see Refs. [tJQII)- It is important to emphasize that the expression 
([2]) is based on a heuristic approximation, which leads to an overestimation of the cross section at high energies ( 
~ 11% at -y/s = 1.3 TeV) in comparison with the more rigorous derivation of the photon spectrum for elastic scattering 
on protons derived in Ref. . For a more detailed comparison among the different photon spectra see Ref. [l2j • As 
a photon stemming from the electromagnetic field of one of the two colliding protons can interact with one photon 
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of the other proton (two-photon process) or can interact directly with the other proton (photon- hadron process), 
both possibilities has been studied in the literature. In principle, the experimental signature of these two processes 
is distinct and it can easily be separated. While in two-photon interactions we expect the presence of two rapidities 
gaps and no hadron breakup, in the inclusive heavy quark photon-hadron production the hadron target we expect 
only one rapidity gap and the dissociation of the hadron. However, as shown in Rcf. |7|, the diffractive heavy quark 
photoproduction, where we also expect the presence of two rapidity gaps in the final state, similarly to two-photon 
interactions, is an important background for two-photon interactions as well as for the dedicated program to search 
evidence of the Higgs and/or new physics in central double diffractive production processes [l3T |. 

In what follows, we briefly review the description of the inclusive and diffractive photoproduction of heavy quarks 
within the dipole picture and the distinct phenomcnological saturation models which will be used in our calculations 
(Section |TTJ) . Our results are compared with the DESY-HERA data 14| and the predictions from the bCGC and 
IP-SAT models for this process are presented for the first time. In Section [IIII the numerical results for the rapidity y 
of the produced states and their total cross sections are shown. A comparison on the order of magnitude of the cross 
sections for the distinct models is performed. Finally in Section Hvl we present our main conclusions. 



II. QCD DIPOLE PICTURE AND SATURATION MODELS 

The photon-hadron interaction at high energy (small x) is usually described in the infinite momentum frame of the 
hadron in terms of the scattering of the photon off a sea quark, which is typically emitted by the small- ir gluons in the 
proton. However, in order to describe inclusive and diffractive interactions and disentangle the small- a; dynamics of 
the hadron wavefunction, it is more adequate to consider the photon-hadron scattering in the dipole frame, in which 
most of the energy is carried by the hadron, while the photon has just enough energy to dissociate into a quark- 
antiquark pair before the scattering. In this representation the probing projectile fluctuates into a quark-antiquark 
pair (a dipole) with transverse separation r long after the interaction, which then scatters off the target [15j . The 
main motivation to use this color dipole approach, is that it gives a simple unified picture of inclusive and diffractive 
processes. In particular, in this approach the inclusive heavy quark photoproduction cross section [jp — > QQX] reads 
as, 

0tot (lP -> QQX) = 2 J d 2 b J d 2 r J dz%(r,z)tf(x,r,b)%(r,z) . (3) 

Furthermore, the diffractive cross section for the process jp — > QQp is given by 

<J? ot {lP^QQp) = j d 2 b J d 2 r J dz^{r,z)N 2 {x,r,b)^> 1 {r,z) . (4) 

In the Eqs. (J3j> and ((4]) the variable r defines the relative transverse separation of the pair (dipole), z (1 — z) is the 
longitudinal momentum fractions of the quark (antiquark) and the function 5 , 7 (r,z) is the light-cone wavefunction 
for transversely polarized photons, which depends in our case of the charge {cq) and mass {tuq) of the heavy quark 
and is given by 

|* 7 (r,z)| 2 =^e 2 Q {[z 2 + (l-z) 2 }m 2 Q K 2 ( mQ r)+ m 2 Q K 2 (m Q r).} . (5) 

The function J\f(x,r,b) is the forward dipole-target scattering amplitude for a dipole with size r and impact 
parameter b which encodes all the information about the hadronic scattering, and thus about the non-linear and 
quantum effects in the hadron wave function (see e.g. (3j). It can be obtained by solving the BK (JIMWLK) 
evolution equation in the rapidity Y = ln(l/x). Many groups have studied the numerical solution of the BK equation, 
but several improvements are still necessary before using the solution in the calculation of observables. In particular, 
one needs to include the next-to-leading order corrections into the evolution equation and perform a global analysis 
of all small x data. It is a program in progress (for recent results see [HI, G3])- ^ n the meantime it is necessary to use 
phenomenological models for Af which capture the most essential properties of the solution. 

During the last years an intense activity in the area resulted in sophisticated models for the dipole-proton scattering 
amp litude, which have strong theoretical constraints and which are able to describe the HERA and/or RHIC data 

on m [n eh ei m, m ee m ei ei . in what follows we will use three distinct phenomcnological saturation models 

based on the Color Glass Condensate which describe quite well the more recent HERA data: the IIM [2l[, the bCGC 
[H E| and the IP-SAT model [H El El- In the IIM model [H[ the scattering amplitude Af(x, r, b) was constructed 
to smoothly interpol between the limiting behaviors analytically under control: the solution of the BFKL equation 
for small dipole sizes, r <C 1/Q s &t{x), and the Levin- Tuchin law [13| for larger ones, r ^> l/Q S at(x). Moreover, the 
authors have assumed that the impact parameter dependence can be factorized: AT(x, r, b) = Af(x, r)S(b). A fit to the 
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structure function i<2(x, Q 2 ) was performed in the kinematical range of interest, showing that it is not very sensitive 
to the details of the interpolation (For details see e.g @). The predictions of this model for several observables were 
studied in Refs. [1, H3, [111, [H)]. Recently the IIM model was improved by the inclusion of the impact parameter 
dependence in the scattering amplitude, with the resulting model being usually denoted bCGC. The parameters of 
this model were fitted to describe the current HERA data in Ref. [24| • Following [22| we have that the dipole-proton 
scattering amplitude is given by: 



M hGGC (x,r,b) = 



■A/n 



1 — exp 



—A In"* (B r ( 



rQ s < 2 
rQ s > 2 



(6) 



with Y = ln(l/x) and k = x"(ls)/x'(ls), where \ is the LO BFKL characteristic function. The coefficients A and 
B are determined uniquely from the condition that Af(x, r) and its derivative with respect to rQ s are continuous at 
rQ a = 2. In this model, the proton saturation scale Q s now depends on the impact parameter: 



Q s = Q s {x,b) = ( — 
\ x 



exp 



b 2 



2-Bcgc / 



(7) 



The parameter Bcgc was adjusted to give a good description of the ^-dependence of exclusive J/ip photoproduction. 
Moreover the factors A/"o and 7 S were taken to be free. In this way a very good description of F2 data was obtained. 
The parameter set which is going to be used here is the one presented in the second line of Table II of (24[: <y s = 0.46, 
Bcgc = 7.5 GeV -2 , A/o = 0.558, xo = 1.84 x 10 -6 and A = 0.119. Furthermore, we will use in our calculations the 
scattering amplitude scattering proposed in Ref. [22j, denoted IP-SAT, which is given by 



AA IP - SAT (x,r,b) 



1 — exp 



2 N, 



■r 2 a s {n 2 )xg{x,p 2 )T{b) 



(8) 



where the scale \x is related to the dipole size r by p = 4/r 

scale up to /i 2 using LO DGLAP evolution without quarks assuming that the initial gluon density is given by 



/i 2 , and the gluon density is evolved from a 
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The values of the parameters /Iq, A g and X g are determined from a fit to F 2 data. 



Moreover, it is assumed that the proton shape function T(b) has a Gaussian form, T(b) = \/{2itBg) exp[— (b 2 /2Bq)], 
with Bq being a free parameter which is fixed by the fit to the differential cross sections for exclusive vector meson 
production. The parameter set used in our calculations is the one presented in the first line of Table III of [22I ]: 
Mo = 1-17 GeV 2 , A g = 2.55, 0.020 and B G = 4 GeV~ 2 . 

As discussed in Refs. [33|, |34| the expression ((5J) for the forward scattering amplitude can be obtained to leading 
logarithmic accuracy in the classical effective theory of the Color Class Condensate formalism. Moreover, it is 
applicable when the leading logarithms in Q 2 dominate the leading logarithms in 1/x, with the small r limit being 
described by the linear DGLAP evolution at small- a;. In contrast, the bCGC model for J\f, Eq. [HI captures the 
basic properties of the quantum evolution in the CGC formalism, describing both the bremsstrahlung limit of linear 
small- ie evolution (BFKL equation) as well nonlinear renormalization group at high parton densities (very small-x). 
Consequently, the IP-SAT model can be considered a phenomenological model for the classical limit of the CGC, 
while the bCGC for the quantum limit. It is important to emphasize that both models provide excellent fits to a wide 
range of HERA data for x < 0.01. Therefore, the study of observables which are strongly dependent on Af is very 
important to constrain the underlying QCD dynamics at high energies. In what follows we consider these two models 
as input in our calculations of the inclusive and diffractive heavy quark photoproduction in jp collisions at HERA 
and coherent pp interactions at Tevatron and LHC energies. For comparison we also consider the IIM model plj . 

Having presented the phenomenological models which will be used in our calculations, in Fig. [T] we compare the 
numerical results for the inclusive heavy quark photoproduction with the experimental DESY-HERA data fl4| ■ In 
all calculations we have use the same quark masses m c = 1.4 GeV and mj, = 4.5 GeV. We quote Ref. Q for a 
comparison of the experimental data with another theoretical approaches. In order to describe the threshold region, 
W — ► 2toq, we have multiplied the cross sections by a factor (1 — a;) 7 , following studies presented in Ref. [35J. In 
charm case (left panel), the IIM and bCGC predictions are almost identical in all kinematical range. In contrast, 
the predictions for these two models for bottom production differ at small energy. The IP-SAT predictions are 
approximately a factor 2 larger than the IIM and bCGC predictions. One have that the IIM and bCGC models 
underestimate the experimental data for charm production at high energies, producing a reasonable description of the 
region near threshold (low energies: W < 20 GeV). In contrast, the IP-SAT model describes the high energy region 
but overestimate the low energy regime. For the bottom case, the value of x which determines the magnitude of the 
saturation scale is not sufficiently small, which implies that the cross section is dominated by the linear regime of the 
scattering amplitude. The three models give a reasonable description of the scarce experimental data. Unfortunately, 
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cc (incl.) 
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3662 nb 


7542 nb 




cc (diff.) 
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161 nb 


532 nb 




bb (incl.) 


51 nb 


51 nb 


158 nb 




bb (diff.) 


0.32 nb 


0.52 nb 


3 nb 



TABLE I: The integrated cross section for the inclusive and diffractive photoproduction of heavy quarks in pp(p) collisions at 
Tevatron and LHC energies. 

the current precision and statistics of the experimental measurements of the photoproduction cross section are either 
low to formulate definitive conclusions about the robustness of the different saturation models presented here. More 
precise measurements could be pose stringent constraints on the energy dependence and overall normalization. Finally, 
it should be noticed that the present calculations concern only the direct photon contribution to the cross section, 
whereas the resolved component has been neglected. In some extent the results from the saturation models presented 
here let some room for this contribution. Details on its calculation and size of its contribution can be found, for 
instance, in Ref. [36j . 

In Fig. [5] we present our predictions for the diffractive photoproduction of heavy quarks. In comparison with 
the inclusive case, the diffractive cross sections are approximately a factor 30 smaller. The main aspect is that the 
difference between the saturations models is enlarged, which is directly associated to the quadratic dependence of the 
cross section on the scattering amplitude. It implies that the experimental study of these observables can be useful 
to determine the QCD dynamics at high energies. 

III. RESULTS 

In what follows, we will compute the rapidity distribution and total cross sections for the inclusive and diffractive 
photoproduction of heavy quarks from proton-proton collisions at high energies. The phenomenological models shortly 
reviewed in the previous section serve as input for the numerical calculations using Eq. |T]) for the energies of the 
current and future pp and pp accelerators. Namely, one considers the Tevatron value \/Snn = 1-96 TeV for its pp 
running and for the planned LHC pp running one takes the design energy V Snn = 14 TeV. 

The distribution on rapidity y of the produced open heavy quark state can be directly computed from Eq. |T]), 
by using its relation with the photon energy u, i.e. y cx ln(w/mg). A reflection around y — takes into account 
the interchanging between the proton's photon emitter and the proton target. Explicitly, the rapidity distribution is 
written down as, 

da [p + p^p + QQ + Y)] JV 7 (w) 

Ty = U ~dT 'tp^** M ' (9) 

where Y is a hadronic final state X resulting of the proton fragmentation in the inclusive case and Y — p for diffractive 
production. 

The resulting rapidity distributions for inclusive and diffractive heavy quark photoproduction coming out of the 
distinct phenomenological models considered in previous section are depicted in Figs. ((3][6|) at Tevatron and LHC 
energies. For the inclusive case (Figs. [3] and [5]) the IIM and bCGC predictions are very similar, as expected from 
the analyzes at photon level in the previous section. In contrast, these predictions are distinct in the diffractive case, 
with the bCGC prediction being larger than IIM one at mid-rapidity. On the other hand, the IP-SAT prediction 
is larger than these predictions by a factor 2 (3) in the charm (bottom) case. We can consider the IIM and bCGC 
predictions as a lower bound for the coherent production of heavy quarks at Tevatron and LHC. Our results indicate 
that the experimental study of the inclusive heavy quark photoproduction can be very useful to discriminate between 
the classical and quantum versions of the CGC formalism. It also is true in the diffractive case, where the different 
models can be discriminated more easily. 

Let us now compute the integrated cross section considering the distinct phenomenological models. The results are 
presented in Table HI for the inclusive and diffractive charm and bottom pair production at Tevatron and LHC. The 
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IP-SAT model gives the largest rates among the approaches studied, followed by the bCGC and IIM models with 
almost identical predictions, as a clear trend from the distribution on rapidity. In the inclusive case, the values are 
either large at Tevatron and LHC, going from some units of nanobarns at Tevatron to microbarns at LHC. Therefore, 
these reactions can have high rates at the LHC kinematical regime. On the other hand, the cross sections for diffractive 
production are approximately two orders of magnitude smaller than the inclusive case, but due the clear experimental 
signature of this process (two rapidity gaps), its experimental analyzes still is feasible. 

In comparison with our previous results for the inclusive production of heavy quarks Q we have that our predictions 
using the modern phenomenological IIM and bCGC models are similar. However, the IP-SAT prediction is factor 
of w 2 larger. In the diffractive case, our predictions are larger by a factor > 2 than those presented in [7|, where 
we have used the GBW model [l8[ as input in our calculations. This behavior is directly associated to the different 
energy dependence predicted by the models for the linear regime. Furthermore, in comparison with the predictions 
for the heavy quark hadroproduction (See e.g. [13]), photoproduction cross sections are smaller (< 1%). However, 
the experimental separation between these two mechanics is feasible due to the presence of one rapidity gap in the 
photoproduction process. 

Lets now calculate the production rates for charm and bottom production in coherent interactions. At Tevatron, 
assuming the design luminosity /Novation = 2 x 10 32 cm s , we have for inclusive production of charm 2 — 4 x 10 2 
and for bottom 2 — 6 events/second. In the diffractive case, we predict 7 - 22 (8 — 60 x 10 -3 ) events/second for charm 
(bottom) production. At LHC, where £lhc = 10 34 cm _2 s~ 1 ), we predict for inclusive (diffractive) charm production 
38 - 75 x 10 3 (16 - 52 x 10 2 ) and for bottom 5 - 15 x 10 2 (3 - 30) events/second. Notice the large rate for bottom at 
LHC. 

Finally, lets discuss the experimental separation of the inclusive and diffractive photoproduction of heavy quarks. 
As emphasized in Refs. [f| 0, E3], although the inclusive photoproduction cross section to be a small fraction of 
the hadronic cross section, the separation of this channel is feasible if we impose the presence of a rapidity gap in 
the final state. It occurs due to the proton which is the photon emitter remains intact in the process. We expect 
that a cut in the transverse momentum of the pair could eliminate most part of the contribution associated to 
the hadroproduction of heavy quarks. Moreover, in comparison with the hadroproduction of heavy quarks, the event 
multiplicity for photoproduction interactions is lower, which implies that it may be used as a separation factor between 
these processes. In the case of diffractive photoproduction of heavy quarks we expect the presence of two rapidity 
gaps in the final state, similarly to two-photon or Pomeron-Pomeron interactions. Consequently, it is important 
to determine the magnitude of this cross section in order to estimate the background for these other channels. In 
particular, the central exclusive diffraction (CED) process characterized by the production of a final state via fusion 
of two Pomerons has being intensely studied as an alternative process to search evidence of the Higgs and/or new 
physics [l3j], with the main background being the exclusive bb production. In Ref. [38|, the double diffractive (DD) 
heavy quark production is studied using the diffractive factorization theorem, including absorption corrections. The 
magnitude of the cross section is the following: for Tevatron one has a^P ~ 4.6 /zb and aPP ~ 0.1 /ib, whereas for 
the LHC one obtains app ~ 18 /xb and aPP ~ 0.5 fib. It is expected that emerging protons from CED and DD 
processes have a much larger transverse momentum than those resulting from diffractive photoproduction processes. 
Consequently, in principle it is possible to introduce a selection criteria to separate these two processes. However, this 
subject deserves more detailed studies. 

IV. CONCLUSIONS 

In summary, we have computed the cross sections for inclusive and diffractive photoproduction of heavy quarks in 
pp and pp collisions at Tevatron and LHC energies, respectively. This has been performed using modern phenomeno- 
logical models based on the Color Glass Condensate formalism, which describe quite well the inclusive and exclusive 
observables measured in ep collisions at HERA. The obtained values are shown to be sizeable at Tevatron and are 
increasingly larger at LHC. The feasibility of detection of these reactions is encouraging, since their experimental 
signature should be suitably clear. Furthermore, they enable to constrain the underlying QCD dynamics at high 
energies, which is fundamental to predict the observables which will be measured in central hadron-hadron collisions 
at LHC. 
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FIG. 1: The total photoproduction cross section for charm (left panel) and bottom (right panel). The experimental measure- 
ments are from DESY-HERA. 




100 
W [GeV] 



10"' 
1000 10 



100 
W [GeV] 



1000 



FIG. 2: The energy dependence of the diffractive cross section for the charm (left panel) and bottom (right panel) production 
predicted by the distinct phenomenological models. 
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FIG. 3: The rapidity distribution for the inclusive charm (left panel) and bottom (right panel) photoproduction on pp reactions 
at Tevatron energy ^Snn = 1-96 TeV. Different curves correspond to distinct phenomenological models. 




FIG. 4: The rapidity distribution for the diffractive charm (left panel) and bottom (right panel) photoproduction on pp 
reactions at Tevatron energy V Snn = 1.96 TeV. Different curves correspond to distinct phenomenological models. 




FIG. 5: The rapidity distribution for the inclusive charm (left panel) and bottom (right panel) photoproduction on pp reactions 
at LHC energy \/Snn = 14 TeV. Different curves correspond to distinct phenomenological models. 
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FIG. 6: The rapidity distribution for the diffractive charm (left panel) and bottom (right panel) photoproduction on pp 
reactions at LHC energy V Snn = 14 TeV. Different curves correspond to distinct phenomenological models. 



